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We study the growth rate and saturation level of the turbulent dynamo in magnetohydrodynamical 
simulations of turbulence, driven with solenoidal (divergence-free) or compressive (curl-free) forcing. 
For models with Mach numbers ranging from 0.02 to 20, we find significantly different magnetic field 
geometries, amplification rates, and saturation levels, decreasing strongly at the transition from 
subsonic to supersonic flows, due to the development of shocks. Both extreme types of turbulent 
forcing drive the dynamo, but solenoidal forcing is more efficient, because it produces more vorticity. 



The turbulent dynamo is the most important process 
to amplify a small initial magnetic field [1] . The growth of 
the magnetic field is exponential, which leads to dynam- 
ically significant magnetic energies on short time scales. 
Dynamo action ranges from the Earth and the Sun [2], 
over the interstellar medium to whole galaxies [3]. Al- 
though the physical conditions (e.g., the different com- 
pressibility of the plasmas) and flow geometries are ex- 
tremely different across these objects, dynamo action has 
been confirmed in all of them. For instance, in the Earth 
and the Sun, the dynamo is driven by subsonic flows. In 
contrast, interstellar clouds and galaxies are dominated 
by highly supersonic, compressible turbulence. 

The main objective of this Letter is to investigate 
fundamental properties of turbulent dynamo amplifica- 
tion of magnetic fields by making systematic numeri- 
cal experiments, in which we can control the compress- 
ibility of the plasma by varying the Mach number and 
the energy injection mechanism (forcing) of the turbu- 
lence. We consider flows with Mach numbers ranging 
from M. = 0.02 to 20, covering a much larger range than 
in any previous study. Haugen et al. [4] provided crit- 
ical Reynolds numbers for dynamo action, but did not 
investigate growth rates or saturation levels, and studied 
only 0.1 < M < 2.6. The energy released by, e.g., super- 
nova explosions, however, drives interstellar and galac- 
tic turbulence with Mach numbers up to 100 [5]. Thus, 
much higher Mach numbers have to be investigated. It is 
furthermore tempting to associate such supernova blast 
waves with compressive forcing of turbulence [6-8]. Mee 
& Brandenburg [6] concluded that it is very hard to ex- 
cite the turbulent dynamo with such curl-free forcing, 
because vorticity is not directly injected. In this Letter, 
we show that the turbulent dynamo is driven by curl-free 
injection mechanisms, and quantify the amplification as a 
function of compressibility of the plasma. This is the first 
study-to the best of our knowledge-addressing the Mach 
number and forcing dependence of the turbulent dynamo 
in detail. The main questions addressed are: How does 
the turbulent dynamo depend on the Mach number of 
the flow? What are the growth rates and saturation lev- 
els in the supersonic and subsonic regimes of turbulence? 
What is the field geometry and amplification mechanism? 



To address these questions, we compute numerical so- 
lutions of the compressible, nonideal, three-dimensional, 
magnetohydrodynamical (MHD) equations with the grid 
code FLASH [9], 

d t p + V • (pu) = 0, 

d t (pu) + V • (pu<g>u - B®B) + Vp* = V • (2vpS) + pF, 
ft£ + V.[(£ + p*)u-(B.u)B] = 

V • [2vpu ■ S + B x (ryV x B)] , 
d t B = V x (u x B) + t?V 2 B, 
V-B = 0, 

(1) 

where p, u, = p + (1/2) |B| 2 , B, and E = pe int + 
(l/2)p|u| 2 + (l/2) |B| 2 denote density, velocity, pressure 
(thermal and magnetic), magnetic field, and total energy 
density (internal, kinetic, and magnetic). Viscous inter- 
actions are included via the traceless rate of strain tensor, 
Sij = (\/2){diUj+djUi) — (l/3)5ijV-u, and controlled by 
the kinematic viscosity, v. We also include physical dif- 
fusion of B, which is controlled by the magnetic diffusiv- 
ity, rj. The MHD equations are closed with a poly tropic 
equation of state, p = c 2 p, such that the gas remains 
isothermal with constant sound speed c s . To drive tur- 
bulence with a given Mach number, we apply the forcing 
term F as a source term in the momentum equation. The 
forcing is modeled with a stochastic Ornstein-Uhlenbeck 
process [8, 10], such that F varies smoothly in space and 
time with an autocorrelation equal to the eddy-turnover 
time, t e d = L/(2A4c s ) at the largest scales, L/2 in the 
periodic simulation domain of size L. M = u rms /c s de- 
notes the root-mean-squared (rms) Mach number, the 
ratio of rms velocity and sound speed. The forcing is 
constructed in Fourier space such that kinetic energy is 
injected at the smallest wave numbers, 1 < |k| L/27T < 3. 
We decompose the force field into its solenoidal and 
compressive parts by applying a projection in Fourier 
space. In index notation, the projection operator reads 

v% (k) = QVtj + (i - Qvl = CSij + (1 - 2C) W|k| 2 , 

where and vfj are the solenoidal and compressive 
projection operators. This projection allows us to con- 
struct a solenoidal (divergence-free) or compressive (curl- 
free) force field by setting £ = 1 (sol) or £ = (comp). 
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For most of the simulations, we set the kinematic vis- 
cosity v and the magnetic diffusivity r] to zero, and thus 
solve the ideal MHD equations. In this case, the dis- 
sipation of kinetic and magnetic energy is due to the 
discretization of the fluid equations. However, we did 
not add any artificial viscosity. Here, we use Riemann 
solvers, which capture shocks also in the absence of arti- 
ficial viscosity. In addition to the ideal MHD simulations, 
however, we also solved the full, nonideal MHD system, 
Eq. 1, for four representative models to show that our 
results are physical and robust against changes in the 
numerical scheme. For the ideal MHD simulations, we 
use the positive-definite, split Riemann scheme HLL3R 
[11] in FLASH v2.5, while our nonideal MHD simulations 
were preformed with the unsplit staggered mesh scheme 
in FLASH v4 [12], using a third-order reconstruction, 
constrained transport to maintain V • B = to machine 
precision, and the HLLD Riemann solver [13]. We ran 
simulations with 128 3 , 256 3 , and 512 3 grid cells, showing 
convergence of our results below. 

We start our numerical experiments by setting L = 
1.24 xlO 19 cm, uniform u = 0, p = 1.93 x 10" 21 gem -3 , 
c s = 2 x 10 4 cms- 1 , and B (0,0, B 0z ) with B 0z = 
4.4 x 10 -16 Gauss in z-direction, corresponding to an ex- 
tremely high initial plasma j3 = 2p/B 2 = 10 20 . These 
values are motivated by dynamo studies of primordial 
clouds [14-16], but in the following, we scale all quantities 
to dimensionless units to address fundamental questions 
of magnetic field amplification in compressible plasmas. 

After an initial transient phase that lasts for 2t e d, tur- 
bulence becomes fully developed and the Mach number 
reaches its preset value, fluctuating on a 10% level. Fig- 
ure 1 (top) shows the time evolution of M in all runs. 
Note the drop in M. for the solenoidally driven runs with 
A4 < 1 as soon as they reach saturation. For these runs, 
the magnetic field has increased to a dynamically signif- 
icant level, causing M to drop at late times, due to the 
back-reaction of B onto the flow. In contrast, in all super- 
sonic runs and in all runs with compressive forcing, the 
magnetic field has little dynamical impact on the turbu- 
lent flow. Although the Mach numbers are not strongly 
affected in those cases, the fragmentation behavior of the 
gas might still change [17], emphasizing the importance 
of magnetic fields. Figure 1 (bottom) shows that the 
magnetic energy grows exponentially over at least 10 or- 
ders of magnitude in each model and reaches saturation 
at different levels (discussed in detail below). Note that 
the nonideal MHD models at different resolution are al- 
most indistinguishable from the ideal MHD models. 

Figure 2 shows that the high Mach number runs 
are dominated by shocks. Compressive forcing yields 
stronger density enhancements for similar Mach num- 
bers [18]. The magnetic field occupies large volume frac- 
tions with rather unfolded, straight field lines in the 
compressively driven cases, while solenoidal forcing pro- 
duces more space-filling, tangled field configurations, sug- 
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FIG. 1. (Color online) Mach number, M (top) and magnetic 
energy _c/m / E m o (bottom) as a function of eddy-turnover time, 
ted, for all runs with solenoidal (sol) and compressive (comp) 
forcing. The Mach number is indicated in the legend. We 
also add nonideal MHD models with M ~ 0.4, 2.5 for sol. and 
comp. forcing, evolved on 256 3 , and 512 3 grid cells. However, 
these models are hardly distinguishable from the correspond- 
ing ideal MHD models, because they are very similar. Thin 
dashed lines show fits in the exponential growth phase. 



TABLE I. Parameters in Eq (3) for the fits in Fig. 3. 
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gesting that the dynamo is more efficiently excited with 
solenoidal forcing. This is quantitatively shown in Fig. 3 
(top and middle panels), where we plot the growth rates, 
T, in the relation E m = £? m oexp(rt), and the saturation 
level, (^ m /^k) sat w ^h the magnetic and kinetic energies 
E m and E^ as a function of Mach number for all mod- 
els. Both T and (E m /E^) sat depend strongly on M and 
on the turbulent forcing. Solenoidal forcing gives growth 
rates and saturation levels that are always higher than in 
compressive forcing, as indicated by the different field ge- 
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FIG. 2. (Color online) Three-dimensional renderings of the gas density on a logarithmic scale for 0.5 < p/po < 50 (from white 
to dark blue), and magnetic field lines (orange) for solenoidal forcing at M = 0.1 (a) and M = 10 (c), and compressive forcing 
at M — 0.1 (b) and M — 10 (d). The stretch-twist-fold mechanism of the dynamo [1] is evident in all models, but operates 
with different efficiency due to the varying compressibility, flow structure, and formation of shocks in the supersonic plasmas. 



ometries shown in Fig. 2. Both Y and (E m / 'E^) sat change 
significantly at the transition from subsonic to supersonic 
turbulence. We conclude that the formation of shocks at 
A4 « 1 is responsible for destroying some of the coher- 
ent vortical motions necessary to drive the dynamo [4]. 
However, as M. is increased further, vorticity generation 
in oblique, colliding shocks [19, 20] starts to dominate 
over the destruction. The very small growth rates of the 
subsonic, compressively driven models is due to the fact 
that hardly any vorticity is excited. To quantify this, we 
plot the solenoidal ratio, i.e., the specific kinetic energy in 
solenoidal modes of the turbulent velocity field, divided 
by the total specific kinetic energy, \ = E so \/E tot in 
Fig. 3 (bottom), which shows a strong drop of solenoidal 
energy for low-Mach, compressively driven turbulence. 
In the absence of the baroclinic term, (1/ ' p 2 )V p x Vp, 
the only way to generate vorticity, w = Vxu, with com- 
pressive (curl-free) forcing is via viscous interactions in 
the vorticity equation [6]: 

ftw = Vx(uxa;) + vV 2 u + 2zA7 x (<SV In p) . (2) 

The second term on the right hand side of the last equa- 
tion is diffusive. However, even with zero initial vorticity, 
the last term generates vorticity via viscous interactions 
in the presence of logarithmic density gradients. The 
small seeds of vorticity generated this way are exponen- 
tially amplified by the non-linear term, V x (ux w), in 
analogy to the induction equation for the magnetic field, 
if the Reynolds numbers are high enough [21]. For very 
low Mach numbers, however, density gradients start to 
vanish, thus explaining the steep drop of dynamo growth 
in compressively driven turbulence at low Mach number. 
Analytic estimates [22] suggest that Y oc M 3 in com- 
pressively driven, acoustic turbulence [23], indicated as 
dotted line in Fig. 3. The solid lines are fits with an 
empirical model function, 

^ = ( po S! +P5 ) xP6 - (3) 



The fit parameters are given in Table I. We emphasize 
that the fits do not necessarily reflect the true asymp- 
totic behavior of Y and (E m /Ek) Siit . The subsonic, 
solenoidally driven models show very high saturation lev- 
els, (E m /E^) sat ~ 40-60%, explaining the strong back 
reaction of the field, causing M. to drop in the satura- 
tion regime (see fig. 1, [24]). For the growth rate, we 
fixed pq such that T oc M 1 ^ 3 for M. ^> 1, in good agree- 
ment with our models up to M « 20. However, even 
higher M has to be investigated to see, if Y oc M 1 ^ 3 
holds in this limit. We find that Y depends much less 
on M. in the solenoidal forcing case than in the com- 
pressive one. Nevertheless, a drop of the growth rate 
at M « 1 is noticeable in both cases. Theories based 
on Kolmogorov's [25] original phenomenology of incom- 
pressible, purely solenoidal turbulence predict no depen- 
dence of T on M.. For instance, Subramanian [26] de- 
rived T = (15/24)Re 1/2 ^ d 1 based on Kolmo gorov-Fokker- 
Planck equations, in the limit of large magnetic Prandtl 
number, Pm = v/r] = Rm/Re ^> 1 with the kinetic 
and magnetic Reynolds numbers Re = Lu rms / (2u) and 
Rm = Lu rms /(2r]). For Pm & 2 [applicable to ideal 
MHD, see 27], and Re « 1500, corresponding to our sim- 
ulations, however, we find slightly smaller growth rates, 
in agreement with analytic considerations [28], and with 
numerical simulations of incompressible turbulence for 
Pm w 1 [29, 30]. Thus, an extension of dynamo theory 
to small Pm is needed. Moreover, extending the theory 
from Kolmogorov to Burgers-type, shock-dominated tur- 
bulence would be an important step forward in develop- 
ing a more generalized theory of turbulent dynamos, po- 
tentially with predictive power for the supersonic regime 
and for compressive turbulent energy injection. 

In summary, we conclude that the growth rate and 
saturation level of the dynamo depend sensitively on the 
Mach number and the energy injection mechanism of 
magnetized turbulence, exhibiting a characteristic drop 
of the growth rate at the transition from subsonic to 
supersonic turbulent flow. Geophysical and astrophys- 
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FIG. 3. (Color online) Growth rate (top), saturation level 
(middle), and solenoidal ratio (bottom) as a function of Mach 
number, for all runs with solenoidal (crosses) and compres- 
sive forcing (diamonds). The solid lines show empirical fits 
with Eq. (3); see Table I. The arrows indicate four models 
(M ~ 0.4, 2.5 for sol. and comp. forcing), using ideal MHD on 
128 3 grid cells (a), nonideal MHD on 256 3 (b), and 512 3 grid 
cells (c), demonstrating convergence for the given magnetic 
Prandtl, Pm^2, and kinematic Reynolds number, Re~1500. 

ical dynamos operate in both, subsonic and supersonic 
plasmas, driven by vastly different injection mechanisms. 
Here we showed that strong magnetic fields are generated 
even in purely compressively (curl-free) driven turbu- 
lence (applicable to e.g., galactic clouds), but solenoidal 
(divergence-free) turbulence drives more efficient dy- 
namos, due to the higher level of vorticity generation 
and the stronger tangling of the magnetic field. 
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